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We report an experimental measurement of a light wavelength at which the ac electric polariz¬ 
ability equals zero for ®^Rb atoms in the F = 2 ground hyperfine state. The experiment uses a 
condensate interferometer both to find this ‘tune-out’ wavelength and to accurately determine the 
light polarization for it. The wavelength lies between the D1 and D2 spectral lines at 790.032388(32) 
nm. The measurement is sensitive to the tensor contribution to the polarizability, which has been 
removed so that the reported value is the zero of the scalar polarizability. The precision is fifty times 
better than previous tune-out wavelength measurements. Our result can be used to determine the 
ratio of matrix elements | ( 5 P 3 / 2 I |d| 1551/2)/(SA/zIMI |55i/2) = 1.99221(3), a 100-fold improvement 

over previous experimental values. New theoretical calculations for the tune-out wavelength and 
matrix element ratio are presented. The results are consistent with the experiment, with uncertainty 
estimates for the theory about an order of magnitude larger than the experimental precision. 

PACS numbers: 03.75.Dg,37.25.-l-k,42.50.Wk 


I. INTRODUCTION 

The energy shift experienced by an atom in an off- 
resonant optical field has found numerous applications 
in atom trapping, manipulation, and measurement. The 
light shift can be characterized by a frequency-dependent 
polarizability, which itself depends in detail on the wave 
function of the electrons in the atom. Accurate mea¬ 
surements of the polarizability can therefore be used to 
test atomic theory calculations, or as phenomenological 
inputs to improve those calculations. Polarizability mea¬ 
surements have a long history of improving our knowl¬ 
edge of atoms in this way mm- 

Precise measurements of the polarizability at optical 
frequencies are technically difficult, because the light 
shift depends also on the optical intensity and it is hard 
to accurately determine the intensity in situ. However, it 
is possible to instead measure a light wavelength at which 
the polarizability equals zero [31]. Since these tune-out 
wavelengths are independent of the intensity, they can be 
accurately measured by various methods |5HS]- 

Tune-out wavelengths can be useful for applications in¬ 
volving species-specific optical manipulation laiiiiiKTo] 
and optical Feshbach resonances [3]. In addition, it was 
recently shown that tune-out wavelengths can be used 
with an atom interferometer for sensitive detection of ro¬ 
tations and accelerations m- Improved knowledge of 
tune-out wavelengths can lead to better performance in 
all these applications. 

In this paper we report measurement of the tune-out 
wavelength for ®^Rb near 790 nm, with an accuracy of 
about 30 fm. This can be compared to the 1.5 to 2 
pm precision of previously reported values for this [5] 
or other tune-out wavelengths laiiiii]. Our result de¬ 
termines the ratio of the D-line dipole matrix elements 
to an accuracy of 15 ppm, about a factor of 100 bet¬ 


ter than previously known [Mill. At our precision the 
measurement is sensitive to many new effects including 
hyperfine interactions [H], QED effects [HI, the Breit 
interaction and the details of the atomic core and 
core-valence interactions [4]. The theoretical tools re¬ 
quired to handle these challenges are closely related to 
those needed for interpreting results such as atomic par¬ 
ity violation and electric dipole measurements in terms 
of fundamental particle properties [18]. Related calcula¬ 
tions are also useful for constraining black-body radiation 
shifts in atomic clocks |19j . Our measurement can thus 
serve as a useful test for theories, or could be taken as a 
phenomenological input value for improved results. 

II. EXPERIMENTAL METHOD 

For an alkali atom in state i, the polarizability can be 
expressed as 

aii^) = T X 2 + Qe + (1) 

h ^ ujfj: - 

where the sum is over all excited states / of the valence 
electron. The transition frequency between i and / is 
ujif and dif = (/|d • e\i) is the dipole matrix element be¬ 
tween i and / for light with polarization vector e. The 
Uc term is the polarizability contribution from the core 
electrons while acv expresses the effect of core-valence in¬ 
teractions [3]. At most frequencies, ac and Ucv are small 
compared to the valence contribution. However, tune- 
out wavelengths occur between pairs of states where the 
valence contributions largely cancel. Figure 1(a) shows 
the tune-out wavelength between the D1 and D2 lines of 
Rb. 

Our measurement uses a Bose condensate atom inter¬ 
ferometer, similar to that previously described in |20) . 
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FIG. 1: Schematic of measurement, (a) Theoretical plot of 
the polarizability a for ®^Rb near the D1 and D2 transitions. 
The polarizability crosses zero at the tune-out wavelength Aq. 
(b) Optical schematic for the experiment. The two Bragg 
laser beams form a standing wave that is used to split and 
recombine a Bose condensate to form an atom interferometer. 
The Stark laser beam illuminates one of the wave packets in 
the interferometer to produce a phase shift, (c) Composite 
image of the atomic wave packets (dark) and the Stark beam 
(white). Here the wave packet centers are 130 /rm apart after 
5 ms of propagation. In (b) and (c), the coordinate axes a:, y 
and 2 are illustrated. 


A condensate of about 10"^ ®^Rb atoms is produced and 
loaded into a weak magnetic trap with harmonic oscil¬ 
lation frequencies of 5.1, 1.1, and 3.2 Hz along the a:, 
y, and z directions, respectively. The trap uses a time¬ 
orbiting potential, with a bias field of 20.0 G rotating in 
the xz plane at 12 kHz frequency. Oscillating magnetic 
gradients provide support against gravity as well as trap 
confinement. 

The atom interferometer is implemented using an off- 
resonant standing-wave laser propagating along the y 
axis, having wave number k. Via Bragg scattering, a 
short pulse from this beam can split the atoms into 
two wave packets traveling with momentum ±2hk HD- 
After 10 ms, the wave packets are reflected using an¬ 
other pulse of the Bragg laser, now adjusted to drive 
the \+2hk) o \—2hk) transition. After 20 ms a second 
reflection pulse is applied, and after another 10 ms, a re¬ 
combination pulse is applied. By using this symmetric 
trajectory, both packets traverse identical paths in the 
trap, which reduces phases shifts and fldelity loss from 
the trapping potential ppj . 

The recombination pulse brings a fraction Nq /N of the 
atoms back to rest in the center of the trap. We obtain 
Nq/N = [1 + V cos(0 -I- (j)r)]/2, where (j) is the phase dif¬ 
ference developed by the atoms during their separation, 
(j}r is the phase shift of the recombination pulse relative 
to the initial splitting pulse, and V « 0.7 is the visibil¬ 
ity. We here set (j)r = tt /2 to maximize the sensitivity to 
(j). We measure Nq /N by allowing the three output wave 
packets to separate for 40 ms and then observing them 
via absorption imaging. 
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FIG. 2: Sample data. The two inset graphs show interfero¬ 
metric measurements of the phase shift (j) induced by the Stark 
beam with power P. The triangles show individual measure¬ 
ments, which are fit to a line to determine the slope. The 
large graph shows how the slope varies as a function of the 
Stark laser wavelength, with the inset graphs corresponding 
to the indicated points. The vertical error bars are the lin¬ 
ear regression errors from the slope flts. The horizontal error 
bars are the standard deviation of several wavelength mea¬ 
surements made over the course of the slope measurement. 
The line in the large graph is another linear flt, and the in¬ 
tercept is taken as our measurement result for the tune-out 
wavelength Aq. Here Aq = 790.03232 nm, with a regression 
error of 50 fm. 


To obtain the polarizability a, we focus another laser 
beam, traveling along z, onto one arm of the interferom¬ 
eter. This Stark beam is applied for 20 ms at the start 
of the interferometer, so that one packet passes through 
it twice. Figure 1(b) shows the orientation of the beams 
involved, and Fig. 1(c) is a composite image of the atoms 
and Stark beam together. 

The energy shift U due to the Stark beam is 


U = -la{£^) 


al 

CqC 


( 2 ) 


where £ is the electric field of the beam, I is the inten¬ 
sity, and c is the speed of light. The brackets denote time 
averaging of the optical field. The light shift induces a 
phase (j) = —{1/h) JUdt proportional to the integrated 
intensity experienced by the atoms. We use an approx¬ 
imately Gaussian beam with waist w « 30 /im. For a 
Stark beam power of P, this yields (j>/{aP) « 66 rad/W 
for a in atomic units. 

The Stark power can be varied from zero to 15 mW us¬ 
ing an acousto-optic modulator. The basic experimental 
procedure is to set the Stark laser to a given wavelength 
A and run the interferometer for different beam powers. 
The resulting phase is fit to a line to determine the slope, 
as shown in the Fig. 2 insets. By performing the exper¬ 
iment at different wavelengths, we plot the slope as a 
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function of A. A second linear fit yields the wavelength 
Ao at which the slope and thus a equals zero. 

III. LIGHT POLARIZATION EFFECTS 

A major complication is that a depends strongly on the 
optical polarization of the Stark beam and the orientation 
of the atomic spin. In general the energy shift can be 
expressed as m 



u = - 


{E^ 


f (0) M ( 1 ) 

-< ’ — Vcos 


+ 


3 cos^ C ~ 1 


3m| - F{F + 1) ( 2 ) 
F{2F-1) “ 


(3) 


where the are irreducible components of the polar¬ 
izability, namely the scalar (i = 0), vector (i = 1) and 
tensor (i = 2) parts. The atom is assumed to be in a par¬ 
ticular hyperfine state iF^mp) relative to the trap mag¬ 
netic field direction b — 'B/B. Here we have F = mp = 2. 
The angle between the Stark beam wave vector and the 
magnetic field is x, so cos x = k-b. Similarly, cos ^ = e-b 
is the projection of the light polarization vector e onto 
the magnetic field. Finally V is the fourth Stokes param¬ 
eter for the light, characterizing the degree of circular 
polarization and expressible as Vcosy = *(e* x e) • 6. 

We are primarily interested in the scalar polarizability 
The tensor contribution is small but measurable, 
and will be discussed below. However, the vector contri¬ 
bution can be quite large. For instance, for <7+ polarized 
light (V = — 1 and x = 0) the vector term completely 
eliminates the tune-out wavelength between the D1 and 
D2 transitions, since the light does not couple our ground 
state to any states in the D1 manifold. To measure the 
tune-out wavelength of the scalar term with the desired 
accuracy, it is necessary to keep |Vcosx| < 10“®. This 
is challenging since it is comparable to the performance 
of the best linear polarizers, and much below the level 
of polarization that can typically be maintained when a 
laser beam passes through a vacuum chamber window. 

We use two methods to control the vector shift. First, 
the rotating bias held of the TOP trap causes cosy in 
Eq. (§ to alternate sign, with a time average close to 
zero. We verihed that the measurement results did not 
depend on the phase of the TOP held at the start of the 
interferometer. 

Second, we linearized the light polarization using the 
interferometer itself. Prior to taking a data set such as in 
Fig. 2, we ran the experiment with the Stark beam pulsed 
on and off synchronously with the TOP held. In this way 
the cosx term could be made close to -1-1 or -I. We ad¬ 
justed the light polarization so that the measured phase 
shifts for those two cases were equal. The polarization 
was established with a calcite polarizer, a zero-order half 
wave plate, and a zero-order quarter wave plate. The 
wave plates could be set to an accuracy of about 0.1°, 
corresponding to V « 2 x 10“^. 


FIG. 3: Effect of tensor polarizability. Data points show 
tune-ont wavelengths Ao, as a function of the angle 0 between 
the linear polarization of the Stark beam and the x axis of 
the trap. Each point is an average of several measurements 
performed as in Fig. 2. For each measurement i the linear fit 
error is combined with the estimated polarization error de¬ 
scribed in the text. The error bars shown are then calculated 
as The solid curve is a sinusoidal fit with a 

variable offset and amplitude. The dashed curve is a fit with 
the amplitude constrained to the expected value. 


After taking the data set, the polarization check was re¬ 
peated and any difference between the cos x = ±1 phases 
was used to estimate the polarization drift that occurred 
during the run. This was converted to a wavelength error 
using an empirical calibration, and the polarization error 
was added in quadrature to the regression error calcu¬ 
lated as in Fig. 2. 

The tensor term in ( ^ gives rise to a dependence on the 
angle of the linear light polarization with respect to the 
trap field, which can be seen in Fig. 3. The polarization 
was adjusted using the half-wave plate in the Stark beam. 
For our geometry, the polarization angle 0 is related to ^ 
in (§ via (cos^ = 0.5 cos^ 0, where the brackets denote 
a time average for the magnetic held. 

Near the tune-out wavelength, a and can be accu¬ 
rately approximated as linear functions {da/dX){X — Aq) 
and {da^^^/dX){X — X^^^) respectively. Here Aq is the mea¬ 
sured value shown in Fig. 3 and A*-®^ is the desired zero 
of the scalar term. The tensor contribution to da/dX is 
negligible, so the two derivatives are nearly equal. If we 
use this and set to zero, we obtain 


in the case of V(cosx) = 0. Fitting to this form, we 
obtain A^^^ = 790.032439(35) nm and {da^^^dX) = 

390(120) fm. This ht is shown as the solid curve in Fig. 

3. 

Alternatively, and da^^^/dX are almost entirely 
due to contributions from the 5P manifold, and can be 
calculated relatively precisely. The derivative term can 














4 


be determined from [15| 


a 


( 0 ) _ 
5P ~ 


10 \dj'\^ uj' 

h^/^^ w'2 - w2 


(-1)^+^' (2F' + 1) 


r 2 1 F' 
\ 1 2 0 


F' 3/2 J' \ 
1/21 2 j 


(5) 


for the F = 2 ground state. Here the sum is over the 
angular momentum quantum numbers of the 5P states, 
w' is the transition frequency to the \ J',F') state, and 
dj' = ( 5 Pj/||d|| 5 «S'i/ 2 ) is the reduced dipole matrix ele¬ 
ment. The dji are known to about 500 ppm precision 
from lifetime and photoassociation measurements | 12 h 

m]. 

Similarly, the tensor term is given by |15j 


( 2 ) _ 
Ogp — 


20 


V 

\dj> 

|V 


Cli'2 _ 

- uj' 

J',F' 



r 

2 1 

F' 

n 

1 2 

2 


(-l)^'( 2 F'-kl) 


F' 3/2 J' 
1/2 1 2 


( 6 ) 


Evaluating the ratio gives /{da/d\) = 538.5(4) fm, 

which is larger than the value determined from our fit by 
about 1 . 3 ( 7 . If we constrain the fit to use the calculated 
value for /{da/dX), we obtain = 790.032388(29) 

nm, about Icr different from the unconstrained result. 
The constrained fit gives a x^/j.o.f. of 1.2, compared 
to 0.5 for the unconstrained fit, both of which are rea¬ 
sonable. Since the calculated value for a^"^^/{da/dX) is 
expected to be accurate, we report the value obtained 
from the constrained fit. 


IV. ERROR ESTIMATION 

As noted, each run of the experiment yields a statisti¬ 
cal error derived from the linear fits of (j) vs. intensity, and 
a polarization error based on the measured polarization 
drift between the start and end of the run. Each run takes 
several hours, so we are not confident that the polariza¬ 
tion change is linear, or even monotonic, throughout the 
run. We therefore use the full value of the polarization 
drift as an error estimate. The polarization drift is in 
fact the largest error source in the measurement. The 
average polarization drift error is 126 fm, compared to 
the average statistical error of 60 fm. Averaging over 
the 21 measurements used would reduce these values by 
However, both errors vary considerably from run 
to run, so for the analysis we combine the two errors for 
each data point in Fig. 3. The resulting fits have the 
uncertainties cited above. 

Another error contribution is the calibration uncer¬ 
tainty in our wavelength measurement. We used a Bristol 
Instruments model 621A wave meter that displayed dig¬ 
its to 1 fm, with results repeatable to about 10 fm. We 
tested the meter by measuring four known saturated ab¬ 
sorption lines in K, Rb, and Cs. The results indicated 


a calibration correction of -40(5) fm at a wavelength of 
790 nm. This correction was applied to the data reported 
here. The full wavelength calibration was performed both 
at the start and end of data collection, and the two Rb 
lines were checked periodically throughout the experi¬ 
ment. No significant differences were observed. 

A significant source of error is asymmetry in the Stark 
laser spectrum [Si- The laser diode source produces 
broadband ASE light [22]. This could be observed 
through its effect on the spontaneous emission rate of the 
atoms, and indicated a background spectral density near 
an atomic resonance of S' « F x 10“^^ Hz“^, in terms of 
the total Stark power P. This is large enough to shift Aq, 
depending on the spectral distribution. We controlled the 
effect by by spectrally filtering the beam using a diffrac¬ 
tion grating and pinhole. Using a 0.4 nm « 200 GHz 
filter bandwidth, and assuming a 10 % variation of the 
spectrum across that bandwidth, the estimated spectral 
density would produce a shift of about 0.1 fm. How¬ 
ever, it is possible that the spectral density near 790 nm 
is larger than that at the atomic resonance. Such low 
spectral power levels are difficult to measure directly, so 
we quantified the effect by comparing our Aq results ob¬ 
tained with 0.2 nm and 0.4 nm filter bandwidths. About 
half our data was taken in each configuration. No mea¬ 
surable difference was observed, within our 30 fm preci¬ 
sion. The expected error would scale as the bandwidth 
squared, indicating that the error for the smaller band¬ 
width configuration was less than 10 fm. We use this as 
the uncertainty from the effect, though we expect it is an 
overestimate. Asymmetry in the tails of the laser line it¬ 
self could similarly shift the measurement, but this could 
be ruled out at the 1 -fm level using an optical spectrum 
analyzer. 

Uncertainty in the trap magnetic field can affect our 
result by changing the value of (cos^ C) in the tensor term. 
The most significant effect is if the magnitude of the 
bias field varies as it rotates. We were able to place a 
limit of 2% on such variations by measuring the Zeeman 
linewidth of the trapped atoms using rf spectroscopy. In 
the worst case, this would induce a 5 fm shift on the 
value of Ai°i. Other effects are smaller, including distor¬ 
tions from a dc background held of less than 1 G, and 
angular misalignment of less than 3 degrees between the 
Stark beam polarization measurement and the plane of 
the bias held. 

The hyperpolarizability of the atoms characterizes the 
nonlinear Stark effect. We estimate the effect by treating 
the Pi /2 and F 3/2 transitions as two-level systems in the 
rotating wave approximation, and summing the resulting 
energy shifts. At the tune-out wavelength, we obtain a 
net shift 


\di/2\^£ 

32h3A3 


where A is the detuning from the P 1/2 transition and £ 
is the Stark held amplitude. At the maximum intensity 
used, this changes A^^^ by only about 1 fm. 
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TABLE I: Estimated error contributions to The entries 
for statistical error and polarization drift report the average 
errors for each type, divided by the square root of the number 
of measurements. In the analysis, both errors were combined 
at each data point to give the reported combined error in the 


result. 

Source Error (fm) 

Statistical 13 

Polarization drift 28 

Stat and polz. combined 29 

Broadband spectrum 10 

Wavemeter calibration 5 

Trap field variation 5 

dc background field 2 

Hyperpolarizability 1 

Total 32 


The effect of interatomic interactions is negligible, as 
the chemical potential of the condensate is only about 
27 rh X 10 Hz. The Zeeman shift from the trap field, how¬ 
ever, is not small. By summing the contributions of the 
individual Zeeman transitions, we calculate that it shifts 
the measured tune-out wavelength blue by 36 fm, so we 
have added this amount to our reported values to give 
the estimated zero-field result. From rf spectroscopy we 
know the bias field magnitude of 20.0(2) G very accu¬ 
rately, so we estimate the error in this shift to be less 
than 1 fm. 

Our error analysis results are summarized in Table |l] 
We sum the errors in quadrature to give our final reported 
one-sigma uncertainty of 32 fm. 


V. COMPARISON TO THEORY 


One other experimental measurement of this tune-out 
wavelength exists, by Lamporesi et al. who obtained 
790.018(2) nm [B]. Our result is in considerable {7a) dis¬ 
agreement, but those authors did not report any special 
effort to control the light polarization. We expect there¬ 
fore that their result is for the particular combination 
of scalar and vector polarizabilities that was relevant to 
their experiment. 

We can however make a useful comparison to theory. 
We Hrst describe how the theoretical result was obtained. 
In the decomposition of Eq. the core terms Oc and a^c 
are are approximately static and are calculated in the 
random-phase approximation [23]. The valence term for 
the 5S state can be expressed in atomic units as 


ay{uj) 


l^ {k\\d\\5S)\Ek-E,s) 

3 “ {Ek — E^s)'^ — 


( 8 ) 


where k = nPi /2 and nP^/ 2 - Up to n = 12 we eval¬ 
uate discrete terms in this sum using experimental val¬ 
ues for the state energies E. Experimental matrix ele¬ 
ments from Ref. [T2] are used for the 5 S' — 6P transitions 


TABLE IL Breakdown of the contributions to the 5S polar¬ 
izability in Rb at A = 790.02568 nm. Reduced matrix ele¬ 
ments d and polarizability contributions are given in atomic 
units. Experimental matrix elements from Ref. [7] are used 
for the 5S — 6 P transitions; remaining matrix elements are 
from the all-order calculations [7l|23]. Uncertainties are given 
in parenthesis. Experimental energies Ai? are measured from 
the ground state and given in cm“^ |25|. 


Contr. 

AE 

d 

ao 

5Pi/2 

12578.951 

4.2199 

-8233.6 

GPi/2 

23715.081 

0.3235(9) 

0.451(3) 

7Pi/2 

27835.05 

0.115(3) 

0.044(2) 

SiPl/2 

29834.96 

0.060(2) 

0.011(1) 

9A/2 

30958.91 

0.037(3) 

0.004(1) 

lOPi/2 

31653.85 

0.026(2) 

0.002 

llF’l/2 

32113.55 

0.020(1) 

0.001 

12Pi/2 

32433.50 

0.016(1) 

0.001 

(n > 12 )Pi /2 



0.022(22) 

5 P 3/2 

12816.54939 

5.9550 

8222.9 

6 P 3/2 

23792.591 

0.5230(8) 

1.173(4) 

7 P 3/2 

27870.14 

0.202(4) 

0.135(6) 

8 P 3/2 

29853.82 

0.111(3) 

0.037(2) 

9 P 3/2 

30970.19 

0.073(5) 

0.015(2) 

IOP 3/2 

31661.16 

0.053(4) 

0.008(1) 

IIF 3/2 

32118.52 

0.040(3) 

0.004(1) 

I 2 P 3/2 

32437.04 

0.033(2) 

0.003 

(n > 12 )P 3/2 



0.075(75) 

Core -1- VC 



8.709(93) 

Total 



0.001 


while all other matrix elements use the all-order calcu¬ 
lations of [23j. The details of the methods are discussed 
in [24]. While experimental values are available for the 
55” — 5P matrix elements m, the theoretical values are 
estimated to have a more accurate ratio, which is most 
important here. For n > 12, the remaining ‘tail’ contri¬ 
butions are calculated in the Dirac-Hartree-Fock approxi¬ 
mation. The state energies and matrix elements are listed 
in Table [ll] Using these values, the tune-out wavelength 
is predicted to lie at = 790.02568 nm as indicated. 


The uncertainty in the theoretical value is dominated 


by uncertainty in the 5P matrix elements. In Table III 


we compare the matrix elements obtained using various 
approximations [23] . All of the methods are intrinsically 
relativistic. The calculations in Table Hill of the tune-out 
wavelength differ only in the values of these two matrix 
elements, with all other values taken from Table [Tl] 


The most accurate methods are expected to be the four 
all-order calculations SD, SDpT, SDgc, and SDpTgc- We 
take the average of these as the final theoretical values, 
and use them to calculate A^°^ = 790.0261(7)nm. The 
uncertainty is estimated from the spread in the four val¬ 
ues. While the scaling (SDsc, and SDpTsc) technique 
is supposed to account for a class of missing correlation 
effects, the scaling affects only about half of the correla¬ 
tion correction in this transition. Therefore we use the 
full spread of the values as our error estimate to allow for 
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TABLE III: Reduced electric-dipole matrix elements for the 55 — 5Pj transitions m, values of the tune-out wavelength 
the matrix element ratio R. Theoretical methods: DF is the lowest-order Dirac-Hartree-Fock, II and III are second- and 
third-order many-body perturbation theory values, SD and SDpT are ab initio all-order values calculated in the single-double 
approximation and with inclusion of the partial triple contributions, and SD^c, SDpTsc are corresponding scaled all-order 
values. Experimental values are averages of several experimental measurements | 12H14 |. 



DF 

II 

III 

SD 

SD,, 

SDpT 

SDpT,, 

Expt. 

55 — 5Pi/2 

4.8189 

4.5981 

4.1855 

4.2199 

4.2535 

4.2652 

4.2498 

4.233(2) 

55 - 5 P 3/2 

6.8017 

6.4952 

5.9047 

5.955 

6.0031 

6.0196 

5.9976 

5.978(4) 

A(°) (nm) 

790.02603 

790.03155 

790.02380 

790.02568 

790.02636 

790.02632 

790.02607 

790.031(6) 

R 

1.9922 

1.9954 

1.9902 

1.9914 

1.9919 

1.9918 

1.9917 

1.995(3) 


TABLE IV: Reduced electric-dipole matrix elements for the 
5S—5Pj transitions and the corresponding line strength ratio 
R. 




DF 

DF-I-Breit 

DF-tQED 

55- 

5Pi/2 

4.8189 

4.8192 

4.82038 

55- 

5 P 3/2 

6.8017 

6.8023 

6.80384 

R 

1.9922 

1.9923 

1.9923 


the effects missed by scaling. We note that this uncer¬ 
tainty estimation is approximate since we are attempting 
to account for unknown correlation effects due to triple, 
quadrupole, and higher excitations. 

The estimated uncertainty in a from all of the non- 
5P contributions is about 0.12 au. Via the derivative 
dX/da = —397 fm/au, this leads to a wavelength error of 
50 fm, about ten times smaller than the uncertainty from 
the 5P levels. The net value of the non-5P contributions 
does give a significant shift of -4.2 pm, mainly from the 
core polarizability. 

The wavelength value determined above does not in¬ 
clude the effects of hyperfine structure. This can be in¬ 
corporated using Eq. ([^ for the 5P levels, using the theo¬ 
retical estimate for the dipole matrix elements. The effect 
of hyperfine structure from all other levels is negligible. 
This yields = 790.0312(7) nm, in reasonable agree¬ 
ment with the experimental value of 790.032388(32) nm. 
The values differ by 1.7cr, with the theoretical uncertainty 
about twenty times larger than that of the experiment. 

As noted, the 5P matrix elements themselves con¬ 
tribute primarily through their ratio 


carry out the DF calculation with the Breit interaction 
included on the same footing with the Coulomb interac¬ 
tion (see, for example. Ref. [H])- The resulting values are 
listed in the column labeled “DF-|-Breit.” Then, we carry 
out the DF calculation with the inclusion of the QED 
model potential, constructed as described in Ref. [16]. 
The Breit interaction is excluded in this calculation to 
separate the two effects. We find that both Breit and 
QED corrections are five times smaller than our uncer¬ 
tainty in the correlation contribution to the ratio. How¬ 
ever, we include the shifts in our estimate R = 1.9919(5). 

An experimental determination of the matrix element 
ratio requires some theoretical input |^. The scalar po¬ 
larizability can be expressed 

= A + \di^ 2 \'^ {Ki /2 + K 3 / 2 R) (10) 

where A includes Oc, cicv, and contributions from valence 
states above 5P. Using the values from Table |Tl| gives 
A = 10.70(12) au. The experimental value for di /2 is 

4.233(2) au [T^HTl] . The coefficients Kji = 

can be obtained from Eq. ([^ and our result for 
Setting = 0 and solving for R yields 1.99221(3). This 
differs from our theory result by O.Gcr, and is about twenty 
times more accurate. Both values are consistent with 
the ratio of the previous experimental matrix elements, 
R = 1.995(3). 


VI. CONCLUSIONS 


R = 


|(5P3/2lMI|5^)P 

|(5Pi/2||d||55)P 


(9) 


This is useful, because the theoretical accuracy of the 
ratio is better than that of the individual matrix elements 
since a large fraction of the correlation corrections cancel. 
This can be seen in the calculations in Table |III[ Using 
the same error estimation procedure as above, we obtain 
a ratio R = 1.9917(5). 

None of the matrix element values in Table ITlll include 
Breit or QED corrections. We evaluate the importance of 
these effects in the lowest-order DF approximation and 


summarize the resulting values in Table IV First, we 


Our measurement of the ®^Rb 790 nm tune-out wave¬ 
length illustrates that tune-out wavelength spectroscopy 
can provide high precision information about atomic ma¬ 
trix elements. As one immediate application, our mea¬ 
surement of the matrix element ratio provides a moder¬ 
ate improvement to the absolute values of the 5Pi/2 and 
5 P 3/2 matrix elements. Each of these elements has been 
determined with about 0.1% precision in three previous 
investigations [HHH]. Our 15 ppm determination of the 
ratio allows all six measurements to be combined, reduc¬ 
ing the total estimated error in each element by about a 
factor of V 2 , as seen in Fig. The resulting best values 
are di /2 = 4.2339(16) and ^ 3/2 = 5.9760(23). Precise 
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FIG. 4: Matrix element values. The points show various de¬ 
terminations of the dj = ( 5 Pj||d|| 5 S'i/ 2 ) matrix elements, for 
J = 1/2 (squares) and J = 3/2 (circles). The first three points 
are measurements by Volz et al. [12] , Simsarian et al. [SI and 
Gutterres et al. M- The fourth point is the error-weighted 
average from the three groups. The fifth point (hollow) is 
the error-weighted average of all six measurements, with the 
constraint ^ 3 / 2/^172 = R ~ 1.99219 obtained in the present 
work. 

knowledge of R may permit yet further improvements 
using the technique of Ref. [26| . 

Our results have several important conclusions in re¬ 
gards to the atomic theory calculations. First, the good 
agreement between the measured and calculated values 
of provides confirmation of the theoretical accuracy. 
Prior to our measurement, the theory result was about 
five times more accurate than the best experimental esti¬ 
mate, making the theoretical prediction difficult to check. 
In particular, our result validates the procedure used to 
estimate the theoretical error, since the error accurately 
reflects the disagreement with experiment. This type of 
error validation is valuable since theoretical error esti¬ 
mates are both challenging and important to obtain. 

Second, we demonstrate that the ratio of matrix ele¬ 


ments can be a useful measure of the accuracy of theo¬ 
retical approaches to include electron correlations. Note 
that the second-order and third-order values in Table uni 
are outside of the theory uncertainty estimate, and dis¬ 
agree significantly with the experimental result. These 
methods are thus confirmed to be less accurate than the 
all-order techniques. 

Third, the accuracy of the experimental ratio value is 
sufficient to test the Breit and QED effects if a more ac¬ 
curate treatment of correlations is carried out. It may be 
possible to achieve this in the full triple coupled-cluster 
approach used to treat Cs parity violation |37|. If suc¬ 
cessful, this would help support the theoretical methods 
and thus clarify the parity violation results dani]. 

The method we have demonstrated can readily be ap¬ 
plied to other tune-out wavelengths in Rb, which we hope 
to pursue in future work. We hope in this way that the 
Rb atom can be established as a well-known reference 
atom for testing theoretical techniques. The significant 
advance in experimental precision should provide a useful 
benchmark for some time to come. 

More generally our method can be applied to any Bose- 
condensed atomic species, which includes many species 
used in precision measurement applications. We hope 
that the improved knowledge of matrix elements made 
possible will prove valuable. 
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